Bioresource Technology 102 (2011) 2047-2052 



Thermal degradation mechanisms of wood under inert and oxidative 
environments using DAEM methods 
D.K. Shen a b , S. Gu a *, Baosheng Jin b , M.X. Fang c 

a Energy Technology Research Group, School of Engineering Science, University of Southampton, Highfield, Southampton S017 1BJ, United Kingdom 

b School of Energy and Environment, Southeast University, Nanjing 210096, China 

c The State Key Laboratory of Clean Energy Utilization, Zhejiang University, Hangzhou 310027, China 


ARTICLE 


INFO 


ABSTRACT 


Article history: 

Received 28 June 2010 

Received in revised form 18 September 

2010 

Accepted 20 September 2010 
Available online 26 September 2010 


Keywords 

Pyrolysis 


Kinetic model 
DAEM 


The pyrolytic behavior of wood is investigated under inert and oxidative conditions. The TGA experiment 
is given a temperature variation from 323 to 1173 K by setting the heating rate between 5 and 40 K/min. 
The results of DTG curves show that the hemicellulose shoulder peak for birch is more visible under inert 
atmosphere due to the higher content of reactive xylan-based hemicellulose (mannan-based for pine). 
When oxygen presents, thermal reactivity of biomass (especially the cellulose) is greatly enhanced due 
to the acceleration of mass loss in the first stage, and complex reactions occur simultaneously in the sec¬ 
ond stage when char and lignin oxidize. A new kinetic model is employed for biomass pyrolysis, namely 
the distributed activation energy model (DAEM). Under inert atmosphere, the distributed activation 
energy for the two species is found to be increased from 180 to 220 kj/mol at the solid conversion of 
10-85% with the high correlation coefficient. Under oxidative atmosphere, the distributed activation 
energy is about 175-235 kj/mol at the solid conversion of 10-65% and 300-770 kj/mol at the solid con¬ 
version of 70-95% with the low correlation coefficient (below 0.90). Comparatively, the activation energy 
obtained from established global kinetic model is correspondingly lower than that from DAEM under 
both inert and oxidative environments, giving relatively higher correlation coefficient (more than 
0.96). The results imply that the DAEM is not suitable for oxidative pyrolysis of biomass (especially for 
the second mass loss stage in air), but it could represent the intrinsic mechanism of thermal decompo¬ 
sition of wood under nitrogen better than global kinetic model when it is applicable. 

© 2010 Elsevier Ltd. All rights reserved. 


1. Introduction 

Bioenergy has attracted more attention than ever before glob¬ 
ally due to its merits of being renewable, widely distributed and 
carbon neutral. Thermochemical processes can have highly effi¬ 
cient conversion of biomass to gaseous, liquid and solid products. 
Industrial scale implementation of these technologies requires 
in-depth understanding of the reaction kinetics for pyrolysis, gasi¬ 
fication and combustion. Many studies have been reported for the 
thermal and kinetic behavior of biomass under inert conditions 
(Koufopanos et al„ 1989; Blasi, 2008) and oxidative conditions 
(Branca and Blasi, 2004; Safi et al„ 2004; Shen et al., 2009). 
However direct comparisons of the characteristics under both 
conditions are not adequately covered, which leaves the scope 
for further study. The kinetic models for biomass pyrolysis (typi¬ 
cally under inert atmosphere) can be categorized as the pseudo- 
mechanistic model with one or more reactions and the distributed 
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activation energy model (DAEM) (Capart et al., 2004; Cai and Liu, 
2008) which is more accurate and recent approach (Sonobe and 
Worasuwannarak, 2008). The DAEM is examined in this study. 

The DAEM is originally developed by Vand (1943). Since then it 
has been applied to analyze the complex reactions in thermal deg¬ 
radation of coal and activated carbon (Pitt, 1962; Solomon et al„ 
1988). Miura employs the DAEM for different types of coals and 
successfully describes a complex reaction which even has a num¬ 
ber of parallel first-order reactions and each occurs at different rate 
(Miura, 1995). Later Miura improves the method to estimate the 
distributed activation energy (/(E)) and the corresponding 
distribution of frequency factor (A(£)) without the assumption of 
functional form of/(£) and A(£) (Miura and Maki, 1998). Recently, 
the nth-order DAEM is proposed considering the influence of tem¬ 
perature on the frequency factor (Cai and Liu, 2008), which gives 
an excellent fit to the pyrolytic mass loss curves of several types 
of biomass without the assumption that all reactions share the 
same frequency factor (A). The latest study includes the validation 
of the DAEM for biomass and estimation of the distributed activa¬ 
tion energies of 62-169 kj/mol at the conversion of 1-84% for corn 
stalk (Li et al„ 2009). To date, there is no reported work for using 


1 All rights reserved. 












2048 


D.K. Shen et al./Bioresource Technology 102 (2011) 2047-2052 


Nomenclature 


A pre-exponential/frequency factor (1/s) 

1 the first stage of thermal decomposition of wood in air 

E activation energy (J/mol) 

2 the second stage of thermal decomposition of wood 

R Universal gas constant (8.3145 J/mol I<) 

in air 

t time (s) 

oo final or infinite 

T temperature (K) 

V mass of the volatiles (g) 

s solid 

X solid mass (g) 

Greek letters 

oc mass fraction at time t or temperature T 

Subscripts 

0 initial 

fi heating rate (K/s) 


the DAEM to describe the oxidative pyrolysis of biomass. This 
study is to fill such knowledge gap and gain better understanding 
of the pyrolytic characteristics of biomass. Thermo-gravimetric 
analysis (TGA) of two wood species is employed to have compara¬ 
tive evaluation of their thermal performances under both inert and 
oxidative environments. The experimental data at different heating 
rates is correlated to obtain the distribution of activation energies 
and coefficients corresponding to different solid conversions. The 
advantages and applications of this method are further estimated 
by comparing with the established global kinetic model in terms 
of the obtained activation energy and correlation coefficient. 


2. Methods 

2.1. Materials 

The properties and the main chemical compositions of the two 
wood species (namely birch and pine) are shown in Table 1. The 
Van Soest method is employed to analyze the main components 
in the wood. The results show that birch consists of more hemi- 
cellulose and less lignin than pine. According to the information 
in the literature (Prins et al„ 2006), birch and pine are categorized 
as deciduous and coniferous wood respectively. Wood samples in 
the same category have found to exhibit very similar pyrolytic 
behavior during the thermal conversion process (Branca and 
Blasi, 2004), while the differences between the thermal decompo¬ 
sition of these two species will be thoroughly discussed in this 
study. 


Properties of wood t 


in the experiments. 


Proximate analysis (%, dry basis) 
Moisture 4.16 

Volatiles 79.04 

Fixed carbon 16.56 

Ash 0.24 

Ultimate analysis (%. dry basis) 
Carbon 50.35 

Hydrogen 6.33 

Nitrogen 0.17 

Oxygen 38.74 

Sulfur 0.01 

High heating value (J/g) 20,843 

Constituent analysis (%) 

Extractives 
Hemicellulose 
Cellulose 
Lignin 

Inorganic salts 


8.7 

17.70 

47.80 

25.50 

0.30 




3.82 

83.50 

12.40 

0.28 


49.53 

6.26 

0.19 

39.85 

0.07 

20,595 


7.59 

28.97 

53.95 

9.43 

0.06 


2.2. Experimental methods 

TGA experiments of solid fuel are performed by a Mettler 
Toledo TGA/SDTA 895IE thermo analyzer under inert (nitrogen) 
and oxidative (air) atmospheres. The samples are prepared by 
grinding large particles of solid fuel in a mortar and sieving. Tests 
are carried out for the particles smaller than 0.3 mm in size and 
less than 5 mg in mass to avoid the effect of mass and heat transfer 
during the process of thermal decomposition (Blasi et al., 2001) at 
the heating rates of 5, 10, 15, 20, 30 and 40 K/min with a temper¬ 
ature range from 323 to 1173 K. During the experiment the 
furnace of TGA is flushed with 50 ml/min nitrogen or air to main¬ 
tain an inert/oxidative atmosphere during thermal decomposition 
of the particles. 

The mass loss of the sample (TG curve) is recorded by an elec¬ 
tronic balance against the linearly increased furnace temperature 
which is adjustably varied with the intensity of the applied radia¬ 
tive heat flux in the TG analyzer. The temperature of the crucible, 
regarded as the real temperature of the solid, is simultaneously 
measured by the thermocouples fixed on the surface of the crucible 
holder to obtain the SDTA curve, to determine the exothermicity or 
endothermicity of the pyrolytic reactions (Branca and Blasi, 2004; 
Fang et al., 2006). It is worth noting that the mass of crucible for 
holding particles is also varied during the heating process of the 
samples. The real mass loss of the particles (fine TG curve) is ob¬ 
tained by subtracting the varied mass of crucible from the crude 
TG curve which is directly recorded by the electronic balance un¬ 
der the same heating rate and atmosphere. Accordingly the mass 
loss rate of the samples (DTG curve) is calculated by differentiating 
the fine TG curve with the temperature without any smoothing 
correction on the data. The DTG curve without smoothing reflects 
the real process of thermal decomposition of wood, although it 
could be the cause for the divergence or non-linear performance 
when estimating kinetic parameter or correlating experimental 
data in the analytical methods (Liu, 2002; Chen, 2006). 

2.3. The mathematical methods 
2.3.1. The DAEM method 

The distributed activation energy model (DAEM) has been em¬ 
ployed successfully for describing coal pyrolysis (Vand, 1943; 
Miura, 1995; Miura and Maki, 1998). Most recent attempts are to 
apply this model to thermal decomposition of biomass and its 
components (Ferdous et al„ 2002; Cai and Liu, 2008; Li et al„ 
2009). In this study, it is assumed that the whole thermal conver¬ 
sion process of wood under both inert and oxidative atmospheres 
is composed of a set of irreversible single (first-order) reactions 
occurring successively. The model is expressed as: 

1 ^ V = 0{EJ)f{E)dE (1) 
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where: 

&(E, T) = exp J* e~ E/RT dT^ (2) 

and 

V = X 0 -X; V oo =X 0 -X oo (3) 

/(£) is the normalized distribution curve of the activation energy, 
representing the variation of the activation energies of those first- 
order irreversible reactions. Through a numerical estimation by 
Miura (1995), the <t>(E, T) function varies steeply with £ at a given 
temperature T for the heating rate /(. A step function U at an activa¬ 
tion energy E s for <2>(E, T) is postulated to give an approximation for 
Eq. (2) as: 

&(£, T) — U(E- E s ) (4) 

And the Eq. (1) could be simplified to: 

\-~=J^m dE ( 5 ) 


environments. For the same V/V„, at different heating rates (/}) cor¬ 
responding to the occurrence of ith reaction, Eq. (11) develops a lin¬ 
ear relationship between {PIT 2 ) and (J/T) with the slope of (-E/R). 
Consequently, the activation energy £ for ith reaction can be deter¬ 
mined by the slope of the linear-fitting curve from the correlated 
experimental data, while the pre-exponential factor A is estimated 
from the intercept. 

2.3.2. The global kinetic model 

The global kinetic model for presenting thermal decomposition 
of biomass is also categorized as the separate-stage model, where 
the notable mass loss region is represented by a global reaction 
(Bilbao et al., 1997; Colomba, 1998; Fang et al., 2006; Shen et al„ 
2009). With regard to the notable mass loss stage(s), the global 
ldnetic scheme for thermal decomposition of wood under inert 
atmosphere is described as by: 

A(Solid) -> B(Char) + C(Gas) 

while the oxidative decomposition of wood is represented by two 
separate-stage reactions in reference to the work by Liu (2002) as: 


This condition is found to hold approximately when E s was chosen 
at <P(E, T) = 0.58 for many combinations of A and/(E) from the pre¬ 
liminary examinations, giving the expression for determining A cor¬ 
responding to activation energy £: 


0.545/i£ = eE/Rr 
ART 2 


(6) 


together with the following approximate equation for: 


<1>(E, T) = exp 



(7) 


This treatment of approximation is given mathematically at the 
temperature T where the ith reaction occurs for a constant heating 
rate p by Miura and Maki (1998): 

(8, 

It shows that the overall reaction rate could be represented approx¬ 
imately by the rate of the ith reaction at the temperature T while 
only ith reaction is occurring, where V f and V ioo , respectively, are 
the amount of volatiles evolved and the effective volatile for the 
ith reaction at the temperature T. Consequently, the whole reaction 
system of biomass decomposition could be approximated by a set of 
single (first-order) reactions, occurring at different temperatures 
corresponding to the successive mass losses or solid conversions 
for a fixed heating rate. 

Eq. (8) could be integrated for a constant heating rate as 
(Freeman, 1963): 


A(Solid) B(Char) + Cl (Gas) 

(Thefirststage) 

B(Char) -> C2(Gas) + D(Ash) 

(Thesecondstage) 

The mass loss fraction a for the stage i: 

is defined as: 

Xio-Xf 

ai = x^x4 

(12) 

The above pyrolytic reactions under inert and oxidative atmo¬ 
spheres are assumed to be governed by first-order Arrhenius Law 


(Zaror and Pyle, 1986; Fang et al., 2006): 


1F-WW 

(13) 

where: 


k = Aexp(-EfRT) 

(14) 

/(*) = !-* 

(15) 

Eq. (13) could be transformed to: 


da k 

m / 

(16) 

An integration function of Eq. (16) is shown as below: 



(17) 

The right integral term of Eq. (17) is approximated by Caost-Red- 
fern method (Caots and Redfern, 1964) to be: 


1 - Vi/V ioo = exp [~jJ t exp {-Ei/RT^j = exp (~ A j^r e 


The Eq. (9) could be rewritten as: 

In Q0 = In 00 - In [- ln(l - V f /V*,)] - ~ ~ (10) 


Since 1 - Vi/Vioo = <P{E, T) m 0.58 was set for deriving the Eq. (6), Eq. 
(10) could be simplified by (Miura and Maid, 1998) 


In 4 = In + 0.6075 - „ ' 


( 11 ) 


Eq. (11) makes it easier to use the DAEM method and obtain kinetic 
parameters for the set of the single reactions during the thermal 
decomposition process of wood under both inert and oxidative 


£ exp {-E/RT)dT » ^ [l exp (-E/RT) 

(18) 

Substituting Eqs. (15) and (18) into Eq. (17), it i 
integration: 

is obtained after 


(19) 

Where: 


g(a) = - ln(l - a) 

(20) 


As the term of 2 RT/E can be neglected since it is much smaller than 
1, Eq. (19) could be simplified as: 




_£ 1 
R T 


(21) 
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It is well-known that the activation energy E could be obtained from 
the scope (-F/R) of Arrhenius curve (In[g(a,)/T 2 ] against 1/T), while 
the intercept of the line with y-axis is related to the pre-exponential 
factor A. This method is applicable to obtain the kinetic parameters 
for the global reaction(s), which is well-established in the literature. 


3. Results and discussions 

3.1. The activation energy from DAEM and global kinetic model 

The inert and oxidative decomposition of wood can be 
described by the DAEM. while the reaction system for the overall 
conversion is assumed to be a set of the single reactions 
corresponding to the successive solid conversions. The kinetic 
parameters according to the degree of solid conversions under 
inert or oxidative atmospheres are determined by the procedures 
as follows (Miura and Maki, 1998): 


-13.2 



M -14.1 
-14.4 


0.95 . 



-14.7 

-15.0 


0.10 



■ 15 K/min 
• 20 K/min 
□ 30 K/min 
o 40 K/min 


1.3 1.4 1.5 1.6 1.7 1.8 1.9 

1/T* 1000 


• Measure V/V., vs. T relationships at different heating rates. 

• Calculate the values of (yS/T 2 ) and ( 1 IT) at the same V/V^ values 
from V/V. x . vs T relationships at different heating rates. 

• Plot ln(/i/T 2 ) vs (1/T) at the same V/V„, value and determine the 
E (from the slope) and A (from the intercept) values from the 
Arrhenius plots at different V)V r% using the relationship in 
Eq. (9). 

The results for birch from the above procedure are illustrated in 
Figs. 1 and 2. The linear and parallel development for different 
V/Voo from 0.1 to 0.85 at various heating rate in Fig. 1 indicates that 
the process of thermal decomposition of birch under nitrogen can 
be described by a set of similar single reactions (Miura and Maki. 
1998). This confirms the previous result that thermal decomposi¬ 
tion of wood under inert atmosphere is mainly determined by 
the degradation of cellulose. Similarly, the Arrhenius plots of birch 
under air are also linear and in parallel for the first stage (the se¬ 
lected V/Voo from 0.1 to 0.65) which is also dominated by the 
decomposition of cellulose. The development is no longer in paral¬ 
lel for the selected V/V^ from 0.7 to 0.95 in Fig. 2, reflecting a set of 
different chemical reactions occurring consecutively (may simulta¬ 
neously) during the second stage of oxidative decomposition of 
wood. 

The activation energy (E) and correlation coefficient (J? 2 ) for the 
two species under both conditions are shown in Figs. 3 and 4. The 


Fig. 2. Parameters from the DEAM method for birch under air at different heating 



Mass loss (%) 


Fig. 3. The activation energy and correlation factor from the DAEM method for 
birch and pine under nitrogen at different mass losses. 
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activation energies at different mass losses {V/V^) for both species 
under nitrogen change slightly from 180 to 220 kj/mol (except the 
case of mass loss at 5% with regard to thermal degradation of the 
extractives), while the correlation coefficients are all above 0.96 
(Fig. 3). The values are higher than those obtained by Li et al. 
(62-169 kj/mol at the conversion of 1-84% for corn stalk) 
(Li et al„ 2009). It is plausible that the thermal decomposition of 
wood under inert atmosphere probably undergoes a set of single 
(first-order) reactions, which could be expressed by the DAEM. 
The activation energies vary narrowly from 175 to 235 kj/mol cor¬ 
responding to the mass loss from 10% to 60 % for both species in air 
(Fig. 4). During the mass loss region of 70-95%, the activation ener¬ 
gies change remarkably from 300 to 770 kj/mol for both species. 
Moreover, the correlation coefficients drop below 0.9 and even 
reach 0.5 (Fig. 4). It is estimated that the stage corresponding to 
the mass loss from 70% to 100% for biomass pyrolysis under oxida¬ 
tive atmosphere constitutes of heterogeneous reactions (mainly 
char combustion), which might be enhanced by the heating rate 
(Sondreal et al., 2001; Fang et al., 2006; Shen et al., 2009). Since 
the same mass loss region in the second pyrolysis stage of wood 
in air is possibly not controlled by the same nominal (first-order) 
reaction at different heating rates, leading to the low linearity of 
the correlated experimental data at different heating rate. 


3.2. The activation energy from global kinetic model 

The activation energy (£) for the global reaction of the inert and 
oxidative decomposition of wood at different heating rate is shown 
in Table 2 and 3, together with the selected temperature region(s) 
corresponding to the main mass-loss stage(s). It is found that the 
activation energy for the two species under nitrogen varies in a 
narrow range between 100 and 120 kj/mol as the heating rate in¬ 
creases, giving the high correlation coefficient (more than 0.96). 
The value for birch is a little higher than that for pine at the same 
heating rate, reflecting the thermal reactivity of the species under 
inert atmosphere. These results imply that the global kinetic 
scheme is reasonable to describe the inert decomposition of wood. 
In air, the activation energy of the first stage reaction, correspond¬ 
ing to the devolatilization of biomass (mainly cellulose and hemi- 
cellulose), only varies slightly (from 109 to 133 K/mol) as the 
heating rate increases, which is similar as the results under nitro¬ 
gen. Comparatively, the substantial scatter of activation energies 
(from 161 to 282 K/mol) in the second stage reaction for both 
two species at different heating rates is possibly attributed to the 
modified mechanism for char oxidation. Another finding is also 
evident by Munir et al. (2009), observing the notable increase in 
value of E with a change of atmosphere from nitrogen to air. 

It needs to be noted that the correlation development of exper¬ 
imental data is employed to obtain the activation energy for both 
two kinetic model (DAEM and global kinetic model). When apply- 


Table 2 

The kinetic parameters for pine and birch under nitrogen (global kinetic model). 


Species Heating rate Temperature E (kj/ A (s ') R 2 

(K/min) range (°C) mol) 


Birch 10 261-382 

15 260-402 

20 270-394 

30 280-406 

40 284-408 

Pine 10 270-408 

15 272-425 

20 270-450 

30 278-458 

40 278-456 


118.25 3.32 * 10 9 0.9644 

117.15 2.98 * 10 9 * 0.9743 

120.43 6.68*10° 0.9663 

122.76 1.14*10'° 0.9684 
124.03 1.9E*10'° 0.9694 
120.30 3.56*10° 0.9835 

114.81 1.26*10° 0.9835 

101.50 9.02 *10 7 0.9791 

103.34 1.51 * 10 s 0.9790 

103.99 2.32* 10 s 0.9792 


Table 3 

The kinetic parameters for pine and birch under air (global kinetic model). 


Species Heating rate Temperature £ (kj/ A (s ') R * 2 

(K/min) range (°C) mol) 



233-391 120.04 

391-490 161.74 

243- 393 120.71 

393-494.5 248.39 

246-396 124.04 

396-478 282.45 

250-398 129.24 

398-488 275.21 

264-402 133.41 

402-500 255.77 

237-408 109.30 

408-530 179.14 

244- 393 121.11 

393-489 225.65 

250-400 118.67 

400-482 254.62 

252-402 123.61 

402-504 214.94 

260-410 121.60 

410-518 201.65 


7.22 * 10° 0.9857 

1.11 * 10" 0.9698 

1.46 *10'° 0.9752 

2.51 * 10 21 4 0.9768 

3.25 * 10'° 0.9782 

4.43 * 10 2 ° 0.9873 

1.15*10" 0.9796 

1.37 * 10 2 ° 0.9926 

2.95 * 10" 0.9731 

5.01 * 10' 8 0.9881 

5.06* 10 s 0.9842 

1.24 * 10 12 * 0.9823 

1.22 * 10'° 0.9810 

1.93 * 10 16 * 0.9838 

7.59 * 10° 0.9813 

2.72 * 10 18 0.9821 

2.74 *10'° 0.9846 

3.87 * 10 15 0.9839 

1.83 * 10'° 0.9818 

3.56 * 10 14 0.9753 


ing global kinetic model, the correlation coefficient ( R 2 ), reflecting 
the linearity of the correlation development and also the validity of 
the obtained activation energy for the global reaction(s), is found 
to be higher than 0.96 for the two species under both inert and oxi¬ 
dative atmospheres. For DAEM, the value (R 2 ) is found to be higher 
than 0.96 for the set of first-order reactions at different mass losses 
under nitrogen, but the value is less than 0.9 for the set of reactions 
corresponding to the mass loss from 70% to 95% under air. From 
the view of apparent estimation, the global kinetic model is appli¬ 
cable for describing thermal decomposition of wood under both in¬ 
ert and oxidative atmosphere, while DAEM is not suitable for 
represent the wood pyrolysis in air. 

However, it is obvious that the obtained activation energy and 
correlation coefficient for the nominal reaction(s) from the global 
ldnetic model is substantially determined by the selection of tem¬ 
perature region(s) corresponding to the main mass loss stage(s) 

varied with the heating rate. For DAEM, the activation energies 
for the set of first-order reactions corresponding to the different 
mass losses (solid conversions) during the thermal decomposition 
process are not directly related to the heating rate, which are 
determined more objectively than that from global kinetic model. 
This implies that the DAEM could represent the intrinsic mecha¬ 
nism of wood pyrolysis better than the global kinetic model, when 
it is applicable. 


4. Conclusions 

The mechanisms of the thermal decomposition of wood are 

examined using DAEM methods and global kinetic model under 

both inert and oxidative atmosphere. The global kinetic scheme 

describes biomass decomposition by the nominated reaction(s) 

dominated in the main separated mass-loss stage(s). Because of 

the high correlation coefficient (R 2 > 0.9) extrapolated from the 

experimental data, the global kinetic model is estimated to be 

applicable for describing the thermal decomposition of wood un¬ 

der both inert and oxidative conditions. DAEM represents the 

whole reaction system of biomass decomposition by a set of 

first-order irreversible single reactions, giving the high activation 

energy than that from the global kinetic model. DAEM is found 

not suitable for describing the thermal decomposition of biomass 

under oxidative condition, because of the complexity of the sec¬ 
ond-stage decomposition reactions. It needs to be noted that DAEM 













2052 


D.K. Shen et al./Bioresource Technology 102 (2011) 2047-2052 


could represent the intrinsic mechanism of wood pyrolysis better 
than the global kinetic model, since it reflects the whole pyrolysis 
system in terms of consecutive solid conversion in spite of heating 
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